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This  report  describes  a  method  for  inverting  VLF/LF  reflection 
data  to  obtain  ionospheric  electron  densities  for  conditions  in  which 
the  geomagnetic  field  causes  the  propagation  to  be  anisotropic.  It 
extends  an  earlier  inversion  method ,  also  developed  for  Rome  Air  De¬ 
velopment  Center,  that  was  useful  only  for  conditions  in  which  the 
propagation  could  be  assumed  isotropic.  The  isotropic  inversion  method 
was  detailed  in  a  companion  report:  R.  E.  Warren,  E.  C.  Field,  Jr., 
and  C.  R.  Warber,  Calculation  of  Ionospheric  Conductivity  Profiles  by 
Inverting  VLF/LF  Reflection  Data:  I.  Isotropic  Propagation ,  Rome  Air 
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This  report  develops  an  "anisotropic"  method  for  inverting  iono- 
sounder  data  to  calculate  ionospheric  profiles.  It  uses  conversion 
coefficients  as  well  as  reflection  coefficients  and,  therefore,  extends 
an  earlier  "isotropic"  inversion  method  that  neglected  the  geomagnetic 
field  and  was  restricted  to  disturbed  conditions  where  collisional 
effects  dominate. 

The  anisotropic  and  isotropic  inversion  methods  differ  primarily 
in  the  calculation  of  Frechet  kernels,  which  become  extremely  compli¬ 
cated  when  the  geomagnetic  field  is  included.  That  complexity  greatly 
increases  the  computer  running  time  needed  to  obtain  numerical  solu¬ 
tions,  and  large  computers  are  therefore  needed  for  routine  applica¬ 
tion.  The  only  other  difference  is  that  the  anisotropic  method  finds 
the  electron  density  profile,  whereas  the  isotropic  method  finds  the 
conductivity  profile.  When  using  the  anisotropic  method,  therefore, 
it  is  necessary  to  input  an  assumed  collision  frequency  profile. 

We  demonstrate  the  anisotropic  method  by  inverting  computer¬ 
generated  reflection  and  conversion  coefficients  calculated  for  known 
electron  density  profiles.  These  coefficients  serve  as  artificial 
data  analogous  to  those  that  would  be  available  from  ionosounders.  We 
ran  three  test  cases:  an  SFE-disturbed  ionosphere,  where  collisional 
effects  were  dominant;  an  ambient  daytime  ionosphere,  where  collisional 
effects  were  comparable  to  geomagnetic  effects;  and  an  ambient  night¬ 
time  ionosphere,  where  geomagnetic  effects  were  dominant.  Despite 
approximations  made  to  accommodate  our  small  computer,  we  obtained 
excellent  agreement  between  calculated  and  known  profiles  for  the  SPE- 
disturbed  and  ambient  nighttime  ionospheres,  and  fair  agreement  for 
the  ambient  daytime  ionosphere. 

Those  results,  which  should  further  improve  when  a  larger  computer 
is  used,  indicate  that  the  inversion  works  best  when  either  collisional 
effects  or  geomagnetic  effects  are  dominant,  and  is  least  accurate 
when  they  are  comparable.  The  examples  show  clearly  that  the  aniso¬ 
tropic  method  must  be  used  for  undisturbed  conditions.  However,  for 
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SYMBOLS 


A 

A  -  admittance  matrix 

c  *  velocity  of  electromagnetic  waves  in  free  space 
C  *  cos  9 

!  *  components  of  electric  field  vector  E  of  wave 

l  *  earth's  magnetic  field 

z  *  Z0(Vx’  Vy'  Vz) 

f  *  components  of  magnetic  field  vector  V  of  wave 
i  -  yCT 
k  *  w/c 

l  *  direction  cosine  of  earth's  magnetic  field 
m  »  direction  cosine  of  earth’s  magnetic  field 

A 

M  -  susceptibility  tensor  for  ionosphere 
n  -  direction  cosine  of  earth's  magnetic  field 
N  -  electron  density 

A 

R  ■  reflection  coefficient  matrix 
S  *  sin  9 
U  -  1  -  iZ 

x  *  cartesian  coordinate  (z  is  measured  vertically  upward 
and  x  horizontally  in  plane  of  incidence) 

2  2 

X  *  Ne  /Eq^u  ,  where  e  and  m  represent  electron  charge 
and  mass,  respectively 

y  ■  cartesian  coordinate  (z  is  measured  vertically  upward 
and  x  horizontally  in  plane  of  incidence) 

Y  *  ligeH  /mw,  where  e  and  m  represent  electron  charge 
and  mass,  respectively 
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z  ■  cartesian  coordinate  (z  is  measured  vertically  upward 
and  x  horizontally  in  plane  of  incidence) 

Z  -  v/w 

1/2 

q  *  characteristic  impedance  of  free  space,  (y^/Cg) 
a  »  azimuth  of  propagation  (measured  east  of  north) 

6  *  magnetic  dip  angle,  X  -  Tt/2 
q  *  electric  permittivity  of  free  space 

9  *  angle  between  vertical  and  wave  normal  of  incident  wave 
below  ionosphere 

X  *  angle  geomagnetic  field  makes  with  vertical  (in  northern 
hemisphere,  tt/2  <  X  s  l) 

q  ■  magnetic  permittivity  of  free  space 

V  ■  collision  frequency  of  electrons 

id  «  2tt  x  frequency  of  waves 


I .  INTRODUCTION 


This  report  extends  previous  analyses  to  include  geomagnetic  effects 
in  calculations  of  ionospheric  conductivity  profiles.  Warren,  Field, 
and  Warber  [1981]  developed  a  method  for  inverting  VLF/LF  ionosounder 
data  to  produce  such  profiles  for  conditions  in  which  geomagnetic  effects 
are  slight.  Subsequently,  Field,  Warren,  and  Warber  [1983]  applied  the 
method  to  reflection  coefficient  data  measured  by  ionosounders  during  a 
strong  solar  proton  event  (SPE) ,  when  propagation  could  be  assumed  iso¬ 
tropic.  They  presented  physical  criteria  for  determining  the  altitudes 
at  which  the  inversion  is  reliable  and  those  where  it  breaks  down. 

Because  it  neglects  the  geomagnetic  field,  the  "isotropic"  inversion 
method  is  restricted  to  conditions  in  which  all  important  reflections 
occur  at  lower  altitudes — less  than  70  km — where  collisional  effects 
dominate  geomagnetic  effects.  Because  of  that  restriction,  the  method 
is  best  suited  to  strongly  disturbed  conditions,  when  ionospheric  reflec¬ 
tion  heights  are  lower  than  under  normal  conditions.  It  gives  only  mar¬ 
ginal  accuracy  for  undisturbed  daytime  conditions,  and  is  inapplicable  to 
undisturbed  nighttime  conditions.  Specifically,  the  isotropic  inversion 

is  valid  only  when  the  ionospheric  conversion  coefficients  are  much 

* 

smaller  than  the  reflection  coefficients. 

The  present  report  develops  an  "anisotropic"  inversion  method  that 
uses  conversion  coefficients  as  well  as  reflection  coefficients  and, 
therefore,  accommodates  the  geomagnetic  field.  We  assume  the  reader  to 
be  familiar  with  the  work  of  Warren,  Field,  and  Warber  [1981]  and  Field, 
Warren,  and  Warber  [1983]  and  do  not  repeat  definitions,  derivations,  or 
source  citations  given  by  them.  The  first  of  those  analyses  is  Vol.  I 
of  the  present  report. 

The  anisotropic  and  isotropic  inversion  methods  are  conceptually 
identical.  They  differ  primarily  in  the  calculation  of  Frechet  kernels. 


These  coefficients  are  often  called  converted  and  normal  reflec¬ 
tion  coefficients,  respectively,  in  the  literature. 
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which  become  extremely  complicated  when  the  geomagnetic  field  is  in¬ 
cluded.  That  complexity  greatly  increases  the  computer  running  time 
needed  to  obtain  numerical  solutions,  and  large  computers  are  therefore 
needed  for  routine  application.  The  only  other  difference  is  that  the 
anisotropic  method  finds  the  electron  density  profile,  whereas  the 
isotropic  method  finds  the  conductivity  profile.  When  using  the  aniso¬ 
tropic  method,  therefore,  it  is  necessary  to  input  a  known  collision 
frequency  profile. 

Because  of  the  computational  complexity  of  the  anisotropic  calcula¬ 
tions,  we  demonstrate  the  method  with  only  three  nominal  cases.  Even  for 
those  numerical  examples,  the  limited  capacity  of  Pacific-Sierra  Research 
Corporation's  Perkin-Elmer  (PE)  3230  computer  required  us  to  make  certain 
approximations  that  would  be  unnecessary  on  a  larger  machine.  Also,  the 
examples  assume  a  vertical  electric  dipole  (VED)  transmitter  and,  there¬ 
fore,  omit  the  reflection  and  conversion  coefficients,  jjRj_  and  j_R||,  that 
pertain  to  the  horizontal  component  of  an  electric  dipole  antenna. 

We  derive  the  Frechet  kernels  in  Sec.  II  and  certain  other  mathe¬ 
matical  formulas  in  the  appendixes.  Those  tedious  derivations  can  be 
skipped  by  the  reader  interested  mainly  in  the  method's  application. 
Section  III  discusses  the  numerical  methods  used  and  tests  the  aniso¬ 
tropic  method  by  inverting  reflection  and  conversion  coefficients  cal¬ 
culated  for  three  known  profiles.  These  coefficients  serve  as  artificial 
data  analogous  to  those  that  would  be  available  from  ionosounders.  Sec¬ 
tion  IV  summarizes  our  conclusions. 


II.  FRECHET  KERNELS  FOR  ANISOTROPIC  IONOSPHERE 


This  section  finds  the  anisotropic  form  of  the  Frechet  kernels. 
Once  those  kernels  are  known,  we  treat  the  reflection  coefficient  terms 
l|R|l  and  jjRj^  as  two  data  points  for  each  frequency  and  proceed  with 
the  inversion  as  given  in  Vol.  I — except  that  here  we  assume  the  colli¬ 
sion  frequency  to  be  known  and  use  the  inversion  to  find  the  electron 
density  profile.  Electric  and  magnetic  quantities  throughout  this 
report  are  in  mks  units. 


COORDINATE  SYSTEM 

Figure  1  shows  the  electromagnetic  wave  normal  and  the  earth's 
magnetic  field  in  the  coordinate  system  used  in  this  report.  We  assume 
that  the  earth  is  flat,  that  the  wave  propagates  in  the  x-direction, 
and  that  the  wave  normal  makes  an  angle  0  with  the  vertical. 


SUSCEPTIBILITY  MATRIX 

To  establish  our  notation,  we  start  with  the  Maxwell  equations  for 
anisotropic  media  written  in  matrix  form  [Budden,  1961]: 


de 

dz 


-ikTe  , 


where 


Note:  Tha  normal  to  tha  warn  front  Um 
in  tha  x-z  plan*  and  makaa  an  angt* 
6  with  tha  vartkai:  Caeoad, 
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Alao.K,  ia  parpandicuiar  to  tha  warn 
normal  in  tha  x-i  plana,  and  H1  it 
parpandicuiar  to  tha  wava  normal 
and  parallal  to  tha  y-axia. 


a.  Electromagnetic  propagation 
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Here,  the  M  are  elements  of  the  susceptibility  matrix,  which  is 


u2  -  i2?2 

-inYU  -  AmY2 

imYU  -  AnY2 

>N  Y 

M  - ^ - =- 

inYu  -  AmY2 

U2  -  m^2 
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L 
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REFLECTION  COEFFICIENTS 

We  can  now  write  the  equation  for  the  anisotropic  reflection  co¬ 
efficients,  generalized  from  Budden's  equation  for  admittance: 

2i  HR  ~21  ~22~  ^^11 

Si  ,  s  +  S  R  -  RS  -  RS  R  , 
k  dz 


where  the  S  matrices  are  given  by 
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When  we  have  free  space  as  defined  by 


X 
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Eq.  (3)  reduces  to  a  particularly  simple  form,  as  follows.  Equations 
(1)  and  (2)  demonstrate  that,  when  Eq.  (5)  holds,  the  only  important 
elements  of  T  are 


23 


Thus,  Eqs.  (4)  reduce  to 
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11 


A 


2CI  , 


yielding,  for  Eq.  (3), 


The  solution  of  Gq.  (6)  is 


Rij 


2Cki(z-Zg) 


(7] 


To  calculate  R,  we  start  high  in  the  ionosphere  where  only  upgoing 
waves  are  present.  We  integrate  Eq.  (3)  downward  until  we  reach  an 
altitude  where  Eq.  (5)  is  true.  Then  we  use  Eq.  (7)  for  the  lower  alti¬ 
tudes.  Appendix  A  details  the  numerical  solution  of  Eq.  (3). 


DERIVATION  OF  FRECHET  KERNELS 

The  procedure  for  determining  anisotropic  Frechet  kernels  is  par¬ 
allel  to  the  isotropic  procedure  of  Vol.  I,  except  that  we  are  now 
dealing  with  matrices.  We  let  R*  be  the  solution  of  Eq.  (3)  for  elec¬ 
tron  density  N  +‘5N,  and  R  be  the  solution  for  electron  density  N.  If 

A 

we  first  define  6R  by 

A 

R'  ■  R  +  6R  and  assume  S'  ■  S  +  5X  , 


where  6X  is  the  change  in  the  magnetoionic  parameter  X,  then  it  can  be 
shown  that 
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We  can  write  Eq.  (8)  in  a  shorter  form  as 


which  we  can  write  as 


-r~  6  $•?<$$  +  Q6X 
k  az 


The  Frechet  kernel  is  defined  by 


<5#(z)  =  I  6  (z,  z')5X(z')  dz' 


Taking  the  derivative,  we  get 


x-632 

k  dz 


[z 

-G(z,  z)6X(z)  +  y  dz  *  "^T 


G  (z,  z')6X(z')  •  (15) 


However,  we  can  also  combine  Eq.  (14)  with  Eqi  (13),  obtaining 


2i  d  *  q\ 


-  ?(z)  I  G  (z,  z')6X(z')  dz’  +  Q(z)6X(z ' ) 


?(z)G(z,  z')5X(z')  dz’  +  Q(z)6X(z) 


Thus,  comparing  Eqs.  (15)  and  (16),  we  see  that 


-  -j^G  (Z,  Z)  *  Q  (z) 


^■§7  G(z,  z')  -  ?(z )G(z,  z’)  . 


(18) 


Since  Che  reflection  coefficient  data  are  known  at  the  ground,  the 
Frechet  kernel  we  need  for  the  inversion  is  G( 0,  z).  However,  the  dif¬ 
ferential  equation  (18)  can  be  time-consuming  to  solve  for  a  large  number 
of  altitudes.  For  small  coupling  between  the  modes,  we  can  use  the 
following  self-consistent  method.  Equations  (17)  and  (18)  and  the  defi¬ 
nition  of  y  in  Eq.  (12)  yield 


Gi<*>  z">  "  -  n“I>  -  H 


| jf  dd"  ?n<2">  jQi<2’ 


J  dz"  [V12(z")  G2(z",  z')  +^13(z")  G3(z",  ..)] 


it  f  *U 


(z"'  )  dz*" 


(19a) 


G2(z,  z') 


-  2lexp 


-  &  /  d2" 


?22(z”)  )  Q2(z ' ) 


z ' 

J  dz”  ^21(z»)  GjU",  z’)  +  ^24(z")  G4(z”,  z') J 


x  exp 


f  dz'”l!  , 


(19b) 
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We  first  solve  Eqs.  (19a)  and  (19d)  assuming  that  G ^  and  G ^  are  zero. 
Then,  to  calculate  G 2  and  G^*  we  substitute  the  results  Into  Eqs. 

(19b)  and  (19c).  Those  results  are  in  turn  reapplied  to  Eqs.  (19a) 
and  (19d),  and  G^,  G^  are  recalculated.  The  process  is  repeated 
until  G^,  G  2»  Gj>  G  ^  stop  changing. 

Computation  time  can  be  saved  if  we  use  the  solution  of  the  Frechet 
kernel  in  the  free  space  between  the  earth  and  the  ionosphere.  Recall 
that,  starting  from  Eq.  (5),  we  found  that 


Therefore,  Eqs.  (10)  yield 


as  long  as  z'  <  zBQ^,  where  zB0T  is  the  bottom  of  the  ionosphere 

Assume  z'  >  z„n_.  Then 
dU  1 


we  obtain  for  Eq.  (23) 


Thus,  we  must  solve  Eqs.  (19)  only  within  the  ionosphere.  Elsewhere, 


TEST  CASES  USING  KNOWN  PROFILES 


The  three  nominal  cases  presented  here  apply  the  anisotropic  in¬ 
version  method  to  conditions  under  which  geomagnetic  effects  are, 
respectively,  slight,  moderate,  and  strong.  We  assume  a  VED  transmit¬ 
ting  antenna  and  therefore  require  only  reflection  coefficients  nR  1( 
and  conversion  coefficients  j|Rj^  to  perform  the  inversion.  However, 
and  jRj|  can  easily  be  incorporated  to  accommodate  the  horizontal  compo¬ 
nent  of  nonvertical  antennas. 

When  geomagnetic  effects  are  important,  the  particle  densities  and 
collision  frequencies  appear  explicitly  in  the  Frechet  kernels  and  ad¬ 
mittance  matrix.  We  have  assumed  the  collision  frequencies  to  be  known 
(see  Appendix  B)  and  solved  for  electron  density  profiles.  This  pro¬ 
cedure  is  different  from  that  used  in  the  isotropic  inversion  of  Vol.  I, 
where  the  assumed  absence  of  the  geomagnetic  field  permitted  the  Frechet 
kernels  and  wave  admittance  to  be  expressed  solely  in  terms  of  the 
scalar  conductivity,  which  was  calculated  directly. 

We  begin  the  inverse  solution  of  the  example  cases  by  calculating 
l|R||  and  yR^  for  the  three  electron  density  profiles  assumed,  and  here¬ 
after  referred  to  as  "true.”  These  coefficients  serve  as  artificial 
data  analogous  to  those  that  would  be  available  from  ionosounders.  All 
calculations  were  done  for  an  azimuth  of  -90  deg,  a  geomagnetic  dip 
angle  of  60  deg,  and  a  geomagnetic  field  strength  of  0.54  G.  Those 
parameters  represent  mid-latitude  east-west  propagation  in  the  northern 
hemisphere.  We  further  assume  an  incidence  angle  of  60  deg. 


COMPUTATIONAL  DIFFICULTIES 

The  computer  code  that  we  developed  for  the  anisotropic  inversion 
is  too  large  to  be  easily  run  on  the  PE  3230.  However,  it  should  be 
within  the  capacity  of  the  large  computers  that  are  commonly  available 
at  government  facilities  and  that  will  be  used  for  routine  applications 
of  the  anisotropic  inversion.  Therefore,  the  size  of  the  code  should 
present  no  problem  in  the  future.  Nonetheless,  it  was  convenient  to 
run  our  test  cases  on  the  PE  3230  and,  to  do  so,  we  had  to  make  two 
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approximations  to  reduce  computer  running  time.  We  emphasize  that  those 
approximations,  discussed  below,  will  probably  be  unnecessary  when  the 
code  is  run  on  a  large  computer.  Moreover,  since  the  PE  3230  uses  only 
32-bit  words,  we  believe  that  round-off  errors  occurred.  We  expect 
greater  accuracy  when  the  code  is  run  on  a  64-bit  machine. 

The  computer  implementation  of  the  anisotropic  inversion  is  straight¬ 
forward — essentially  the  same  as  described  in  Vol.  I — except  that  the 
calculation  of  the  Frechet  kernels  is  extremely  time-consuming.  Even 
using  the  self-consistent  calculation  described  in  Sec.  II,  we  found  the 
execution  time  on  the  PE  3230  to  be  about  1  CPU-minute  per  frequency  per 
iteration.  A  typical  anisotropic  inversion  requires  about  5  frequencies 
and  30  iterations  and,  therefore,  about  150  CPU-minutes.  Such  execution 
times  are  impractical,  so  we  made  two  approximations. 

Our  first  approximation  was  to  double  the  spacing  between  the  alti¬ 
tudes  at  which  the  kernels  are  calculated,  and  then  interpolate  to  find 
the  kernels  at  the  intermediate  altitudes.  Second,  and  of  greater  con¬ 
sequence,  we  assume  the  kernels  to  quickly  fall  to  zero  at  altitudes 
above  those  where  the  self-consistent  approximation  is  valid.  Such  an 
approach  avoids  exhaustive  calculation  of  the  kernels  at  altitudes  too 
high  to  have  much  effect  on  the  ground-level  reflection  coefficients. 
Numerical  checks  indicate  that  the  inversion  depends  only  slightly  on 
the  value  of  the  kernels  at  such  high  altitudes.  That  fact,  coupled 
with  the  satisfactory  agreement  between  true  and  calculated  profiles 
found  for  the  examples  discussed  below,  indicates  that  both  approxima¬ 
tions  are  good,  although  probably  unnecessary  for  a  large  computer. 

Their  practical  effect  is  a  fourfold  reduction  in  computation  time. 

EXPONENTIAL  FIT  TO  SPE-DISTURBED  PROFILE 

For  our  first  example,  we  assume  an  exponential  electron  density 
profile  given  by 


N  (z) 
e 


2  x 


in9  (z-45) / 4 
10  e 


e/m" 


(25) 


When  combined  with  the  electron  collision  frequency  profile  in  Appen¬ 
dix  B,  Eq.  (25)  produces  an  exponential  conductivity  profile  that 
approximates  the  polar  ionosphere  between  about  40  and  50  km  at  1508  UT 
on  4  August  1972,  the  peak  period  of  a  strong  SPE  [Reagan  et  al.,  1982]. 
That  conductivity  profile  was  used  in  Vol.  1  as  an  example  for  testing 
the  isotropic  inversion  method. 

Table  1  lists  the  reflection  and  conversion  coefficients  calculated 


for  the  profile  given  by  Eq.  (25).  As  expected  for  a  strong  SPE,  the 
conversion  coefficients  are  much  smaller  than  the  reflection  coefficients 
The  anisotropic  inversion  should  therefore  produce  nearly  the  same  re¬ 
sults  as  the  isotropic  inversion. 


Table  1 

REFLECTION  AND  CONVERSION  COEFFICIENTS  CALCULATED 
FOR  MODEL  OF  SPE-DISTURBED  IONOSPHERE,  EQ.  (25) 


Frequency 

(kHz) 

Reflection 

Coefficient, 

l|lRl|l 

Conversion 

Coefficient, 

1  „>tjj 

1  ||Rll/l||R||l 

5 

0.341 

0.009 

7 

0.307 

0.009 

0.030 

10 

0.267 

0.009 

18 

0.167 

0.007 

0.040 

30 

0.069 

0.003 

0.048 

Figure  2  shows  the  agreement  between  the  true  profile  and  that 
calculated  using  the  anisotropic  inversion.  As  expected,  these  results 
correspond  well  to  those  reported  in  Vol.  I  for  the  isotropic  inversion. 
Specifically,  we  find  close  agreement  between  the  true  and  calculated 


The  electron  density  profile  in  Eq.  (25)  is  physically  unreal¬ 
istic  in  two  respects.  First,  an  SPE  would  not  disturb  the  ionosphere 
at  the  mid-latitude  assumed  here.  Second,  ions  rather  than  electrons 
dominate  the  conductivity  below,  say,  50  km,  so  Eq.  (25)  is  actually 
an  "effective"  electron  density.  Nonetheless,  the  profile  is  suitable 
for  our  present  purpose — namely,  verifying  the  mathematical  accuracy 
of  the  anisotropic  inversion. 
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profiles  at  altitudes  from  about  40  to  50  km,  and  poorer  agreement  out¬ 
side  that  range. 

In  Vol.  I  and  the  subsequent  paper  by  Field,  Warren,  and  Warber 
[1983],  we  identified  physical  criteria  for  determining  the  altitudes 
at  which  VLF/LF  reflection  data  yield  reliable  profiles.  The  most  use¬ 
ful  of  those  criteria  was  based  on  calculation  of  an  upgoing  wave,  if, 
normalized  to  unity  at  the  ground.  The  calculated  profiles  were  most  re¬ 
liable  at  altitudes  where  0.5  <,  ff  <  1,  because  the  ionosphere  at  those 
altitudes  most  strongly  affects  the  ground-level  reflection  coefficients. 
At  low  altitudes,  where  ff  «  1,  the  ionosphere  is  too  rarefied  to  cause 
substantial  reflection.  At  high  altitudes,  where  H  s  0.5,  the  signal  is 
too  weak  and  the  return  too  weak,  even  if  the  reflectivity  is  high. 

Inclusion  of  the  geomagnetic  field  precludes  derivation  of  a  simple 
quantity  like  H.  Therefore,  we  derive  in  Appendix  C  two  complicated 
quantities,  U11  and  Uj,,  that  combine  upgoing  electric  and  magnetic 
fields  and  can  be  used  to  assess  the  range  of  validity  of  the  aniso¬ 
tropic  inversion — much  as  S’  was  used  to  assess  the  range  of  validity  of 
the  isotropic  inversion.  The  quantities  and  are  normalized  to 
1  and  0  at  the  ground,  respectively.  However,  does  not  become  if 
if  the  geomagnetic  field  vanishes,  nor  is  the  application  of  and 
for  the  anisotropic  inversion  so  unequivocal  as  the  application  of  H 
for  the  isotropic  inversion.  One  difficulty  is  that  the  differential 
equations  for  have  a  mathematical  property  known  as  stiffness ,  which 
causes  numerical  inaccuracies — but,  fortunately,  only  at  high  altitudes 
that  exceed  those  where  the  most  important  reflections  occur.  We  use 
and  in  all  three  nominal  cases  examined  in  this  report. 

Figures  3  and  4  plot  and  for  the  calculated  profile  shown 
in  Fig.  2  and,  respectively,  frequencies  of  5  and  30  kHz.  For  both 
frequencies,  the  condition  0.5  <,  <  1  is  satisfied  in  the  altitude 

range  for  which  the  calculated  and  true  profiles  agree,  and  is  not 
satisfied  for  higher  or  lower  altitudes.  Therefore,  this  example  con¬ 
firms  that  provides  a  means  of  assessing  the  range  of  validity  of 
an  inversion  based  on  ionosounder  data  for  which  no  true  profile  is 
available  for  comparison. 


50  6 

Altitude  (km) 


Fig.  3— Normalized  field  strength  versus  altitude  for 
calculated  profile  of  Fig.  2  and  5  kHz 


Field,  Warren,  and  Warber  [1983]  used  the  isotropic  inversion 
method  to  calculate  conductivity  profiles  from  ionosounder  data  mea¬ 
sured  during  a  strong  SPE.  That  method  uses  only  j|R|(  data  as  inputs, 
so  jjRj_  data  must  be  ignored.  The  error  incurred  by  that  procedure 
could  not  be  determined,  but  it  was  argued  that  great  accuracy  should 
be  expected  because  ^R^  <<  at  all  frequencies  used. 

We  can  now  evaluate  the  accuracy  of  the  above  procedure  by  using 
values  of  ||R||  from  Table  1  as  inputs  to  the  isotropic  inversion  method 
and  ignoring  The  results,  shown  in  Fig.  5,  are  reassuring.  In 

fact,  the  agreement  between  the  calculated  and  true  profiles  is  slightly 
better  in  Fig.  5  than  for  the  anisotropic  inversion  plotted  in  Fig.  2. 
That  apparently  puzzling  result  is  easily  explained.  The  slight  in¬ 
crease  in  accuracy  gained  by  retaining  ,,Ri  is  more  than  offset  by  the 


Altitude  (km) 


Fig.  4— Normalized  field  strength  versus  altitude  for 
calculated  profile  of  Fig.  2  and  30  kHz 


loss  of  numerical  accuracy  incurred  by  using  the  more  detailed  aniso- 
* 

tropic  method.  Because  for  this  example  the  anisotropic  inversion 
requires  much  longer  computer  time  without  giving  more  accurate  re¬ 
sults,  we  conclude  that  the  isotropic  method  should  always  be  used  if 
„R1  «  ijRji  at  all  frequencies  for  which  data  are  available. 


This  conclusion  applies  to  results  calculated  on  the  PE  3230 
and  could  be  different  for  calculations  on  larger  computers,  for  which 
round-off  errors  are  smaller  and  approximations  are  not  required. 


Altitude  (km) 


Fig. 


-Calculated  and  true  profiles  for  SPE-disturbed  ionosphere 
Eq.  (25)  and  isotropic  inversion,  using  nR,,  as  input 
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EXPONENTIAL  FIT  TO  AMBIENT  DAYTIME  PROFILE 

Our  second  example  assumes  an  exponential  electron  density  profile 
given  by 


Ne<2) 


4  x  108 


g0. 35(z-70) 


(26) 


which  approximates  ambient  daytime  conditions.  The  calculated  reflec¬ 
tion  and  conversion  coefficients  are  shown  in  Table  2.  Since  the  con¬ 
version  coefficients  are  of  the  same  order  of  magnitude  as  the  reflec¬ 
tion  coefficients,  geomagnetic  effects  should  be  important.  The  large 
ratio  of  (|Rj^  to  j|R |(  at  the  higher  frequencies  is  caused  by  ||R||  being 
small  rather  than  jjRj^  being  large.  That  effect,  which  could  be  due  to 
a  quasi-Brewster's  angle,  could  cause  the  geomagnetic  effect  at  the 
chosen  parameter  values  (for  incidence  angle,  azimuth)  to  be  greater 
than  for  other  daytime  conditions.  The  computational  limitations  of 
the  PE  3230  preclude  calculating  several  cases  to  test  that  possibility. 


Table  2 

REFLECTION  AND  CONVERSION  COEFFICIENTS  CALCULATED  FOR 
MODEL  OF  AMBIENT  DAYTIME  IONOSPHERE,  EQ.  (26) 


Frequency 

(kHz) 

Reflection 

Coefficient, 

1  |lRlll 

Conversion 
Coefficient , 

1  |,8xl 

Ii|8ll/I„»||l 

5 

0.277 

0.181 

0.654 

7 

0.233 

0.200 

0.858 

10 

0.194 

0.211 

1.090 

18 

0.133 

0.194 

1.459 

30 

0.092 

0.132 

1.430 

The  numerical  results  for  ambient  daytime  conditions,  shown  in 
Figs.  6  through  9,  follow  the  same  sequence  as  those  of  Figs.  2  through 
5.  Figure  6  shows  that  the  anisotropic  inversion  gives  good  results  at 
altitudes  between  about  64  and  70  km,  and  poor  results  outside  that 
range.  Overall,  however,  the  anisotropic  results  are  poorer  than  those 
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Fig.  6— Calculated  and  true  profiles  for  ambient  daytime 
ionosphere— Eq.  (26)  and  anisotropic  inversion 
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Fig.  7— Normalized  field  strength  versus  altitude  for 
calculated  profile  of  Fig.  6  and  5  kHz 
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-Normalized  field  strength  versus  altitude  for 
calculated  profile  of  Fig.  6  and  30  kHz 


Iculated  and  true  profiles  for  ambient  daytime 
phere— Eq.  (26)  and  Isotropic  inversion, 
using  „R„  as  input 


obtained  with  the  isotropic  method  and  presented  by  Field,  Warren,  and 
Warber  [1983],  which  showed  close  agreement  between  calculated  and 
true  values  in  the  60  to  70  km  range  for  a  nominal  daytime  profile. 

The  degradation  in  accuracy  is  almost  certainly  caused  by  computational 
errors  associated  with  the  complicated  anisotropic  equations,  but  re¬ 
calculation  on  a  larger  computer  is  needed  to  verify  that  explanation. 

Figures  7  and  8  show  the  normalized  upgoing  fields  and 
As  in  the  previous  example,  the  altitude  range  defined  by  the  condition 
0.5  £  <  1 — roughly  60  to  70  km — corresponds  to  the  64  to  70  km 

range  where  the  anisotropic  inversion  works  best. 

Recall  that  the  better  agreement  obtained  by  Field,  Warren,  and 
Warber  [1983]  with  the  isotropic  method  was  based  on  artificial  data 
generated  by  neglecting  the  geomagnetic  field,  and  does  not  imply  that 
the  isotropic  inversion  should  be  used  instead  of  the  anisotropic  in¬ 
version  for  reflection  coefficients  actually  measured  during  the  day¬ 
time.  Figure  9  confirms  the  somewhat  poorer  accuracy  that  would  be 
achieved  by  using  the  isotropic  method  to  invert  measured  reflection 
coefficients.  Specifically,  we  ignore  the  values  of  shown  in 
Table  2  and  use  only  values  of  j|Rjj.  As  expected,  the  agreement  between 
the  calculated  and  true  profiles  is  not  quite  as  good  as  that  shown  in 
Fig.  6  for  the  anisotropic  inversion. 

AMBIENT  NIGHTTIME  PROFILE 

We  used  the  profile  given  by  Pappert  and  Moler  [1974],  labeled 
true  profile  in  Fig.  10,  to  test  the  anisotropic  inversion  under  am¬ 
bient  nighttime  conditions.  The  reflection  and  conversion  coefficients 
calculated  for  that  profile  are  given  in  Table  3.  Since  ||Rj^  is  large, 
geomagnetic  effects  would  be  expected  to  dominate  collisional  effects. 
Figure  10  shows  that  the  calculated  and  true  profiles  agree  very  well 
at  altitudes  between  73  and  85  km.  In  fact,  the  overall  agreement  is 
much  better  than  that  achieved  for  the  ambient  daytime  model  and  indi¬ 
cates  that  the  anisotropic  inversion  works  best — on  the  PE  3230  com¬ 
puter,  at  least — when  geomagnetic  effects  are  dominant. 


Frequency 

(kHz) 


liiRxUli,Ri,l 


Reflection  Conversion 
Coefficient,  Coefficient, 

IiiriiI  !  iirj_I 


5 

0.322 

0.392 

1.218 

7 

0.267 

0.467 

1.746 

10 

0.251 

0.504 

2.011 

18 

0.335 

0.418 

1.250 

30 

0.275 

0.318 

1.159 

Figures  11  and  12  show  that  the  condition  0.5  ^  U  <  1  works 
well  to  define  the  height  range  over  which  the  anisotropic  inversion 
is  most  reliable.  Further,  the  behavior  of  indicates  the  spurious¬ 
ness  of  the  layer  shown  by  Fig.  10  to  occur  at  around  60  km.  That 
layer  affects  the  fields  only  slightly,  because  collisional  effects 
cause  electrons  at  60  km  altitude  to  be  much  less  mobile  than  electrons 
at  73  to  85  km.  Its  artificiality  is  confirmed  by  Figs.  11  and  12, 
which  show  that  changes  only  slightly  in  the  vicinity  of  the  spur¬ 

ious  layer.  It  is,  of  course,  possible  that  the  layer  would  not  appear 
if  the  inversions  had  been  performed  on  a  larger  computer. 

Figure  13  shows  the  consequences  of  ignoring  ||R_[_  and  using  only 
liRjl  in  the  isotropic  inversion  method.  Such  a  procedure  is  equivalent 
to  using  the  isotropic  inversion  method  in  conjunction  with  ionosounder 
data  measured  at  night.  As  expected,  the  calculated  and  true  profiles 
totally  disagree.  Of  course,  the  anisotropic  method  must  be  used 
for  ambient  nighttime  conditions. 


Altitude  (km) 
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IV.  CONCLUSIONS 


We  tested  the  anisotropic  inversion  method  for  three  cases:  an 
SPE-disturbed  ionosphere,  where  collisional  effects  were  dominant;  an 
ambient  daytime  ionosphere,  where  collisional  effects  were  comparable 
to  geomagnetic  effects;  and  an  ambient  nighttime  ionosphere,  where  geo¬ 
magnetic  effects  were  dominant.  Despite  approximations  made  to  accom¬ 
modate  our  small  computer,  we  obtained  excellent  accuracy  for  the  SPE- 
disturbed  and  ambient  nighttime  ionospheres,  and  fair  accuracy  for  the 
ambient  daytime  ionosphere. 

Those  results,  which  should  further  improve  when  a  larger  computer 
is  used,  indicate  that  the  inversion  works  best  when  either  collisional 
effects  or  geomagnetic  effects  are  dominant,  and  is  least  accurate  when 
they  are  comparable.  The  examples  show  clearly  that  the  anisotropic 
method  must  be  used  for  undisturbed  conditions.  However,  for  disturbed 
conditions,  when  )(R£  <<  ^R^,  it  is  preferable  to  neglect  j(Rj^  and  use 
values  of  j|R|j  as  inputs  to  the  isotropic  method,  which,  for  that  situ¬ 
ation,  gives  quite  accurate  results  with  modest  computer  running  time. 
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Appendix  A 

CALCULATION  OF  FORWARD  PROBLEM 


This  appendix  presents  our  method  of  calculating  the  forward 
problem — which  determines  the  reflection  coefficients,  given  the  elec¬ 
tron  density.  Strictly,  the  forward  problem  is  the  solution  of  Eq.  (3) 
in  Sec.  II.  In  practice,  however,  we  can  save  computer  time  by  solving 
for  the  admittance  rather  than  the  reflection  coefficients,  and  then 
using  [Budden,  1961] 


to  find  R.  Since  we  solve  only  for  A  within  the  ionosphere,  we  avoid 
numerical  instabilities. 

DEFINITION  OF  ADMITTANCE 

We  start  with  Budden's  equation  for  admittance  [Budden,  1961]: 


where 


k  3z 


T21  +  AT11 


l22 


12 


(A.  1) 


(A. 2a) 


-T 


14 


•12 


0 


1 


(A. 2b) 


Here,  the  are  elements  of  Budden  s  T  matrix  [Eq.  (1)].  Substituting 
Eqs.  (A. 2)  into  Eq.  (A.l),  we  have  the  forms  for  the  admittance  in  the 
general  case: 


k  dz  A11  "  T41  “  ^^l  +  A11(T44  "  Tll)  +  A12A21  ’ 


k  dz  A12  "  ”T42  +  A12(-T14A11  +  A22  +  T44)  +  T12A11  ’ 


k  dz  A21  “  T31  +  A21('T14A11  +  A22  “  Tll)  +  T34A11  » 


—  ^  A  *  — T  +A^-TAA  +TA  +TA 

k  dz  22  L32  A22  i14A12A21  *12  21  *34*12 


INITIAL  VALUES  FOR  ADMITTANCE 

We  will  use  our  normal  coordinate  system.  Assume  a  sharply  bounded 
ionosphere  with  a  wave  incident  from  below,  and  assume  we  know  q  (the 
solution  of  the  Booker  quartic  in  the  ionosphere) .  Then  the  Maxwell 
equations  become 


qEy  "  ‘Hx  ’  "qEx  +  SEz  *  "Hy  •  _SEy  *  ‘Hz  * 


and 


qH  -  [(I  +  M)E]  ,  -qH  +  SH  -  [(I  +  M)E]  ,  -SH  -  [(1  +  M)E] 

y  a  a  z  y  y 


By  combining  these  equations,  we  obtain 

r  2  . 


2  2 

0  -q  - 


qs 


0  E  ■  -(I  +  M)E 


”S  E 


(A.  3) 


or,  with  obvious  notation, 

A  /V  A  "► 

(L  +  I  +  M)E  *  0 

Then,  for  a  nontrivial  solution  to  exist. 


■  A  A  A  a 

L  +  I  +  M  -  0  . 


This  leads  to  Booker's  quartic,  which  we  solve  using  the  method  of 
Sheddy  [1968]. 

Let  q  be  a  solution  of  Booker's  quartic  and  define 

D(q)  -  L(q)  +  I  +  M  . 

Now  expanding  Eq.  (A. 3),  we  get 


’l2Ey  +  D13EZ  ’ 

0  , 

(A.  4a 

22Ey  +  B23Ez  ‘ 

0  , 

(A.  4b 

32Ey  +  D33Ez  " 

0  . 

(A.  4c 

>  .  <*  /'  i*  .*  •'  V'  '  ■  ’  ’  *  '  »  *  -  W  •  »  •  ,  "  L  '  k  •  -  V  **  L*  •  ”•  ■  *•»  %  \  _  *■  *1  S  a  ' 


To  find  expressions  for  E  and  E  in  terms  of  E  ,  we  use  Eqs.  (A. 4a) 

x  z  y 

and  (A. 4c).  Then 


E  •  -  n33n12  'n-n'2  E  =  PE  , 
X  D33D11  ~  D13D31  y  7 


D31D12  ”  D11D32 

Ez  “  d  D  -  DD  Ey  E  QEy 
°33ull  u13u31  y  7 


(A.  5a 


(A.  5b 


Also, 


Hx  «  -qEy  , 


Hy  “  qEx  “  SEz  "  (qP  "  SQ)Ey  2  ^y 


(A.  6a 


(A. 6b 


The  definition  of  admittance  [Budden,  1961]  is 


(A.  7) 


where  the  superscripted  (1)  and  (2)  denote  one  solution  of  the  Booker 
quartic.  There  are,  in  general,  four  solutions:  two  represent  upgoing 
waves  (Imq  <  0)  and  two  represent  downgoing  waves  (Imq  >  0) .  At  this 
point,  only  upgoing  waves  are  present. 

If  we  let  *  aE^,  then  we  can  solve  Eq.  (A. 7).  Substituting 

Eqs.  (A. 5)  and  (A. 6)  into  Eq.  (A. 7)  gives 

*Here  we  adopt  the  notation  of  Sheddy  [1968]. 


Equation  (A. 8)  works  only  in  the  anisotropic  case.  In  the  iso 
tropic  case. 


where 

M  »  n2  -  1  . 

The  Maxwell  equations  separate  into  two  sets:  the  transverse 
electric,  i.e.. 


qE  -  -H  ,  -SE  -  -H  ,  -qH  +  SH  «  n2E  ; 
y  x  y  z  x  z  y 

and  the  transverse  magnetic,  i.e., 

-,Ex  +  SEz  -  -Hy  ,  ,Hy  -  n2Ex  ,  -SHy  -  n2Ex  . 

2  2  2 

Both  sets  lead  to  q  -  n  -  S  . 


We  can  take 


Appendix  B 


COLLISION  FREQUENCY 


This  appendix  presents  the  electron/neutral  ion  collision  frequency 
as  a  function  of  height.  We  take  the  values  for  collision  frequency 
from  the  tables  given  by  Pappert  and  Mbler  [1974].  We  use  nighttime 
collision  frequency  data  for  all  cases  because  daytime  and  nighttime 
frequencies  are  the  same  up  to  160  km,  and  we  do  not  expect  to  need 
frequencies  above  160  km  for  daytime  cases. 


Table  B.l 

ELECTRON/NEUTRAL  ION  COLLISION 
FREQUENCY  VERSUS  ALTITUDE 


Collision 


Altitude 

(km) 

Frequency 

,  'lx 
(sec  ) 

250.00 

1.05E  02 

225.00 

3.50E  01 

220.00 

3.00E  01 

210.00 

3.30E  01 

200.00 

4.50E  01 

150.00 

1.60E  03 

120.00 

1.00E  04 

100.00 

3.90E  04 

0.0 

4.30E  11 
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Appendix  C 

CALCULATION  OF  UPGOING  FIELDS 


This  appendix  derives  the  equations  for  the  upgoing  fields.  Since 
nearly  all  quantities  here  are  matrices,  we  drop  the  caret  matrix  nota¬ 
tion.  Budden  [1961]  defines  the  up-  and  downgoing  fields  as,  respectively 


U  -  (I1E  +  I2H)  , 
D  -  i  (I3E  +  I4H)  , 

where 


We  also  know  that  H  *  AE  and  D  *  RU. 

Combining  Eqs.  (C.l)  and  (C.2),  we  see  that 


(C.l) 

(C.2) 


(C.  3) 


But  it  is  easy  to  show  that 


Now  using  Eqs.  (C.l)  and  (C.4)  and  the  fact  that 


V31 


> 


we  obtain 


H  -  i"1!!  +  I~LD  . 


(C. 


Budden  gives  the  derivatives  of  E  and  H  as 


f 


-T11E  *  T12H  • 


(C. 


i  H'  '  TnE  +  T22H  ,  (C. 

where  T  is  as  defined  by  Eqs.  (A. 2).  If  we  differentiate  Eq.  (C.l), 
we  get 


£  U'  -  \  (IjE'  +  I2H’)  . 


(C. 


Putting  Eqs.  (C.6)  and  (C.7)  into  Eq.  (C.8)  and  then  using  Eqs.  (C.4) 


After  some  matrix  manipulation,  we  obtain 


and 


(C.l 


(C.l 


The  matrices  and  S2  are  the  same  as  and  defined  in 

Sec.  II. 

Entering  D  »  RU  into  Eq.  (C.9)  yields 

£  U'  -  \  (sx  +  S2R)U  .  (C.l 

It  is  useful  to  know  the  fields  in  free  space.  Section  II  shows  that 

S^  ■  2CI  and  S2  *  0  . 


Therefore,  Eq.  (C.12)  becomes 
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